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ABSTRACT 



We analyze radial velocity observations of the 12 extremely low-mass <0.25 
M white dwarfs (WDs) in the MMT Hypervelocity Star Survey. Eleven of the 
12 WDs are binaries with orbital periods shorter than 14 hours; the one non- 
variable WD is possibly a pole-on system among our non-kinematically selected 
targets. Our sample is unique: it is complete in a well-defined range of apparent 
magnitude and color. The orbital mass functions imply that the unseen com- 
panions are most likely other WDs, although neutron star companions cannot be 
excluded. Six of the 11 systems with orbital solutions will merge within a Hubble 
time due to the loss of angular momentum through gravitational wave radiation. 
The quickest merger is J0923+3028, a g = 15.7 ELM WD binary with a 1.08 
hr orbital period and a <130 Myr merger time. The chance of a supernova la 
event among our ELM WDs is only 1% - 7%, however. Three binary systems 
(J0755+4906, J1233+1602, and J2119-0018) have extreme mass ratios and will 
most likely form stable mass-transfer AM CVn systems. Two of these objects, 
SDSS J1233+1602 and J2119-0018, are the lowest surface gravity WDs ever 
found; both show Ca II absorption likely from accretion of circumbinary mate- 
rial. We predict that at least one of our WDs is an eclipsing detached double 
WD system, important for constraining helium core WD models. 



*Based on observations obtained at the MMT Observatory, a joint facility of the Smithsonian Institution 
and the University of Arizona. 

** Spitzer Fellow 
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SDSS J092345.60+302805.0, SDSS J105353.89+520031.0, SDSS J123316.20+160204.6, 
SDSS J142200.74+435253.2, SDSS J143948.40+100221.7, SDSS J144801. 13+134232.8, 
SDSS J151225.70+261538.5, SDSS J163026.09+271226.5, SDSS J211921. 96-001825.8) 
— white dwarfs 



INTRODUCTION 



Extremely low mass (ELM) WDs, with masses <0.3 M Q , are the remnants of stars that 
never ignited helium in their cores. The Universe is not old enough to produce ELM WDs by 
single star evolution. Thus ELM WDs must undergo significant mass loss during their evo- 
lution. Althoug h metal-rich red g iants with strong winds may evolve in isolation into single 



-0.4 M WDs (IKilic et a 



200 7c|), producin g ~0.2 M Q ELM WDs most likely requires com- 



pact binary systems (e.g. iMarsh et al.lll995l ). Observational data for ELM WDs is limited, 
however, because of their rarity. For example, among the 9316 WDs identi fied in the Sloan 



Digital Sky Survey (SP SS), fewer than 0.2% have masses below 0.3 M Q ( lEisenstein et al. 
20061 : iKepler et al.ll2007h . 



Kilic et al.l ( 120091 . l2010bl ) have established a radial velocity program to search for com- 



panions around known ELM WDs. Of th e six ELM WDs o bserved to date - including 
J0917+4638, the lowest mass WD known kilic et al.l l2007bh - all six ELM WDs are in 
binaries with <1 day orbital periods. Two more recently identified ELM WDs, NLTT 1174 8 



and J1257+5428 f Kawka fc Vennes 2009: Marsh et al. 



are also in binari es with <1 day orbital periods ( jSteinfadt et al. 



2010 



Kulkarni fc van K erkwiik 



2010 



201 



Kawka et al. 



0) 



2010 



Kilic et al.l l2010al ). Three of t hese eight ELM WDs will merge due t o gravitational wave 
radiation in less than 500 Myr ( iMullally et al.l 120091 ; IKilic et al.ll2010bl ). A larger sample of 
ELM WDs is required to measure the space density, period distribution, and merger rate of 
these systems. 

Here we present 12 ELM WDs w it h < 0. 25 M^ fo und in the Hypervelocity Star (HVS) 
Survey of iBrown et al.l (120051 |2006ajJbJ, l2007al Jbl. l2009al lbl). ten of which are new discoveries. 
Our ELM WD sample is unique because it comes from a complete, non-kinematically-selected 
survey of stars targeted in a well-defined range of apparent magnitude and color. Radial 
velocity follow-up reveals that 11 of the 12 ELM WDs are binaries with <14 hr orbital 
periods. Clearly, the compact binary formation scenario is the best explanation for ELM 
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WDs. 

In Section 2 we describe the survey design and the spectroscopic observations. In Section 
3 we derive the physical parameters of the ELM WDs with stellar atmosphere modeling, and 
present the radial velocity curves for each object. In Section 4 we discuss the nature of the 
binary systems, and conclude in Section 5. 



DATA 



2.1. Survey Design 



The ELM WDs reported here are found in the HVS survey, a radial velocity survey of 
objects with the colors of late-B type stars. All survey targets were selected from the Sloan 
Digital Sky Survey (SPSS) photometric ca talog using de-reddened, uber-calibrated PSF 
magnitudes (lAdelman-McCarthy et al.l 120081 ) . The HVS survey color selection is illustrated 
in Figure [TJ Although the color selection was designed to exclude normal WDs by their 
(u — g) color, the selection fortuitously includes low surface gravity WDs (see Figured]). 

The HVS survey consists of three parts. The first part of the HVS survey is defined by 
15 < g < 17 (dash-dot line in Fig, m and is 100% complete over 7300 deg 2 of the SDSS 
Data Release 5 footprint iBrown et al.l ( 2007al ). The second part of the HVS survey is defined 
by 17 < g < 19.5 (dashed line in Fig. [J) and is 100% c omplete over the same 7300 deg 2 of 
the SDSS Data Release 5 footprint IBrown et al.l ( l2007bl ). The final part extends the earlier 
surveys over 9800 deg 2 of the SDSS Data Releas e 7 footprint and also extends the magnitude 
limit to go = 20.5 (dotted line), as described in lBrown et al.l (j2009al ). The final part of the 
HVS survey is currently 88% complete. 

Our spectroscopy reveals that 15% of the HVS S urvey targets are WDs. T wo of the 
WDs with ~0.2 M are previously published elsewhere ( IKilic et al.ll2007al . l2010bl ). Here, we 
present the remaining ten WDs with mass <0.25 M Q . 



2.2. Spectroscopic Observations 

With the exception of J0923+3028, we obtained all observations at the 6.5m MMT 
telescope using the Blue Channel spectrograph. We operate the spectrograph with the 832 
line mm -1 grating in second order, providing wavelength coverage 3650 A to 4500 A and a 
spectral resolution of 1.0 - 1.2 A, depending on the slit size used. We obtain all observations 
at the parallactic angle, with a comparison lamp exposure paired with every observation. 
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Fig. 1. — Color-color diagram showing the observed distribution of ELM WDs (solid stars) 
and the other WDs (crosses, circles) found in the HVS survey (15<g <17 dot-dash line, 
17<go<19.5 dashed line, 19.5<go<20.5 dotted line) compared with our synthetic photometry 
of the hydrogen atmosphere WD models of D. Koester (solid lines). 
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We obtained spectroscopy for J0923+3028 at the Fr ed Lawrence Whipple Observatory 
1.5m Tillinghast telescope using the FAST spectrograph (IFabricant et al.lll998l ). We operate 
the FAST spectrograph with the 600 line mm -1 grating, providing wavelength coverage 3600 
A to 5500 A and a spectral resolution of 2.3 A. All observations are paired with a comparison 
lamp exposure. 

We process the spectra using in the st andard way. We flux-calibrate using blue 

spectrophotometric standards ( Massey et ali 1988), and we measure radial velocities using 
the cross-correlation package RVSAO (IKurtz fc Minklll998l ). 



2.3. Radial Velocities 

It is important to maximize velocity precision in our tests for variability, and we achieve 
the best precision by cross-correlating the ELM WDs with themselves. Our procedure begins 
by cross-correlating the observed spectra with a high signal-to-noise WD template to obtain 
preliminary velocities. We then shift the spectra to the rest frame, and sum them together to 
create templates for each target. Finally, we cross-correlate the spectra with the appropriate 
template to obtain the final velocities for each target. The average radial velocity uncertainty 
of the measurements is ±18 km s _1 . 

We check our velocities using WD model spectra with atmospheric parameters cus- 
tomized for each target. The results are consistent within 10 km s _1 , which we take as 
our systematic uncertainty. Table [3] in the Appendix presents the full set of radial velocity 
measurements for the 10 newly discovered ELM WDs presented here. 



3. RESULTS 



Our time-series spectroscopy provides for robust determinations of effective temperature 
and surface gravity for each object as well as its binary orbit al parameters. Follow-u p 
spectroscopy of J0917+4638 and J1053+5200 is already published fcilic et al-lboOTTjl boiObh . 
but the other ELM WDs in our sample have not been studied until now. 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Fig. 2.— Model fits (cyan lines) overplotted on the composite observed spectra (black 
lines). The spectral continua provide improved T e g determination except for the poorly 
fluxed-calibrated J2119— 0018, for which we use the continuum-corrected spectrum. 
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Fig. 3. — Surface gravity vs. effective temperature of the observed ELM WDs (red stars) and 
the other WDs (crosses, circles) found in the HVS survey, c ompared with predicted tracks 



for He WDs with 0.16-0.45 Mm (b l ue lines. LPanei et al. 



2007) and CO WDs with 0.5-1.0 M Q 



(green lines, iBergeron et al.lll995t iHolberg fc Bergeron! 120061 ). Our sample of ELM WDs is 
defined by m < 0.25 M . 
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3.1. Stellar Atmosphere Parameters 



We perform stellar atmosphere model fits using synthetic DA WD spectra kindly pro- 
vided by D. Koester. The grid of WD model atmospheres covers effective temperatures from 
6000 K to 30,000 K in steps of 500 K to 2000 K, and surface gravities from logg = 5.0 to 9.0 
in steps of 0.25 dex. The model atmospheres are calculated ass uming local th ermodynamic 
equilibrium and include both convective and radiative transport (jKoesterll2008l ). We perform 
fits on the average composite spectra for each object and present the resulting T e fj and log g 
values in Table [TJ We also perform fits to the individual spectra to derive a robust statistical 
error estimate for each object, also presented in Table [TJ Figure [2] shows the composite 
spectra and our model fits. The surface gravities range 5.1-6.9 dex at effective temperatures 
of 10400-18300 K, confirming that the objects are ELM WDs. 

We fit the flux-calibrated spectral continua to better measure effective temperature. 
The exception is J2119— 0018, the one ELM WD we observed under variable conditions 
and at high airmass, for which we fit the continuum-corrected spectrum. If we fit only the 
continuum-corrected Balmer line profiles for all the ELM WDs, we obtain best-fit solutions 
that differ by 440 ± 190 K in T e s and 0.04 ± 0.04 dex in logg from our published values. 
These differences reflect our systematic error, and demonstrate that our fits to the entire 
flux-calibrated spectra are reasonably accurate. Our fits also agree well with the SDSS 
photometry in all five filters, an additional demonstration that our temperature and surface 
gravity measurements are reliable. 

Figure [3] compares the obs erved ELM WDs and other WDs in the HVS survey with 
the improved iPanei et al.l (120071) track s fseelKilic et al.ll2010bl ) for He-core WDs with masses 
0.16-0.45 M and the Bergeron et al. ( 1995 )^1 tracks for normal CO-core WDs with masses 
0.5-1 M . The gap between the 0.17 M and 0.18 M He-core WD tra cks is linked to th e 
threshold for thermonuclear flashes in the hydrogen shell burning phase ( iPanei et al.l 120071 ) . 
Remarkably, the majority of observed ELM WDs fall in the gap of parameter space between 
the 0.17 M and 0.18 M models. 

The inconsistency between the observations and the He-core WD models make accurate 
mass and luminosity estimates difficult. Fortunately, mass and luminosity change slowly 
over the range of effective temperature and surface gravity sampled by our ELM WDs. We 
conclude that the majority of the ELM WDs have a mass of ~0.17 M and an absolute 
magnitudes of M g ~ 8; more precise estimates are possible for the 0.20 - 0.25 M ELM WDs 
and are summarized in Table [TJ 



3 http://www. astro. umontreal.ca/~bergeron/CoolingModels/ 
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Using these absolute magnitude estimates, the majority of our ELM WDs are found 
at heliocentric distances of 1 < d < 3 kpc. The notable exception is J0923+3028, a bright 
g = 15.7 ELM WD approximately 280 pc distant. All of our ELM WDs are located at high 
Galactic latitudes 30° < |6| < 80° in the SDSS imaging footprint. 

The systemic velocities suggest that two objects are halo WDs: J0818+3536 and J1422+4352 
have syste mic heliocentric ra dial velocities of —201 ± 4 km s _1 and —195 ± 8 km s _1 , re- 



spectively. iKilic et al.l (l2010bl ) also identify J1053+5200 as a halo object based on its proper 
motion and distance estimate. Unfortunately, reliable proper motions are unavailable for our 
fainter g ~ 20 WDs. Despite having relatively large > 1 kpc distances above the Galactic 
plane, the remaining ELM WDs in our sample have systemic radial velocities consistent with 
a disk origin. 



3.2. Orbital Parameters 

Eleven of our twelve ELM WDs exhibit significant radial velocity variation, with peak- 
to-peak velocity amplitudes up to 890 km s~ x . We compute orbital periods for these systems 
by finding the period that minimizes x 2 f° r a circular orbit. Figure H] plots the periodograms 
for the 10 new ELM WDs. A few ELM WDs have multiple period solutions because of 
insufficient coverage, however in all cases the periods are constrained to be <1 day. We 
estimate the period error by conservatively identifying the range of periods with \ 2 < ^Xmim 
where Xmin * s the minimum \ 2 ■ 



We compute best-fit orbital elements using the code of iKenyon fc Garcia ( 119861 1 . which 
weights each velocity measurement by its associated error. The uncertainties in the orbital 
elements are derived from the covariance matrix and \ 2 . To verify these uncertainty esti- 
mates, we perform a Monte Carlo analysis where we replace the measured radial velocity v 
with v + gSv, where Sv is the error in v and g is a Gaussian deviate with zero mean and unit 
variance. For each of 10000 sets of modified radial velocities, we repeat the periodogram 
analysis and derive new orbital elements. We adopt the inter-quartile range in the period 
and orbital elements as the uncertainty. For binaries with multiple period aliases, both ap- 
proaches yield similar uncertainties. When there are several equally plausible periods, the 
Monte Carlo analysis selects all possible periods and derives very large uncertainties. In 
these cases, we adopt errors from the covariance matrix for the lowest \ 2 orbital period. 



- 10 - 




Time [days] 



Fig. 4. — Periodograms for the 10 new ELM WDs. Some objects have multiple period 
aliases because of insufficient coverage, however in all cases the periods are constrained to 
<1 day. J1448+1342 is formally consistent with no velocity variability. 



Table 1. White Dwarf Physical Parameters 



Object 


So 


(u - g)o 


(.9 - r)o 


T e s 


logff 


M s dfceiio Mass 




(mag) 


(mag) 


(mag) 


(K) 




(mag) (kpc) (M ) 



J0755+4906 


20 


095 


± 





023 





763 


± 





086 


-0 


304 


± 





039 


13160 


± 


260 


5 


84 


± 





07 


8.0 


2 


62 





17 


J0818+3536 


20 


484 


± 





030 





766 


± 





152 


-0 


274 


± 





058 


10620 


± 


380 


5 


69 


± 





07 


8.0 


3 


14 





17 


J0917+4638 1 


18 


696 


± 





019 





738 


± 





037 


-0 


328 


± 





024 


11850 


± 


170 


5 


55 


± 





05 


8.0 


1 


38 





17 


J0923+3028 


15 


628 


± 





018 





190 


± 





027 


-0 


401 


± 





021 


18350 


± 


290 


6 


63 


± 





05 


8.4 





28 





23 


J1053+5200 2 


18 


874 


± 





023 





395 


± 





049 


-0 


358 


± 





045 


15180 


± 


600 


6 


55 


± 





09 


8.7 


1 


08 





20 


J1233+1602 


19 


829 


± 





018 





809 


± 





068 


-0 


233 


± 





028 


10920 


± 


160 


5 


12 


± 





07 


8.0 


2 


32 





17 


J1422+4352 


19 


794 


± 





020 





715 


± 





061 


-0 


271 


± 





032 


12690 


± 


130 


5 


91 


± 





07 


8.0 


2 


29 





17 


J1439+1002 


17 


812 


± 





016 





434 


± 





042 


-0 


341 


± 





026 


14340 


± 


240 


6 


20 


± 





07 


8.0 





92 





18 


J1448+1342 


19 


217 


± 





023 





581 


± 





047 


-0 


286 


± 





033 


12580 


± 


230 


6 


91 


± 





07 


9.9 





73 





25 


J1512+2615 


19 


241 


± 





021 





551 


± 





051 


-0 


258 


± 





030 


12130 


± 


210 


6 


62 


± 





07 


9.3 





97 





20 


J1630+2712 


20 


040 


± 





019 





818 


± 





083 


-0 


223 


± 





031 


11200 


± 


350 


5 


95 


± 





07 


8.0 


2 


50 





17 


J2119-0018 


20 


000 


± 





021 





867 


± 





092 


-0 


197 


± 





033 


10360 


± 


230 


5 


36 


± 





07 


8.0 


2 


51 





17 



References. — 1 iKilic et all |2007bh : 2 iKilic et al.l (fcoiObh 
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Table 2. Binary Orbital Parameters 



Object 




p 

(days) 




K 

(km s -1 ) 


^systemic 

(km s- 1 ) 


Spec. Conjunction 
(days + 2450000) 


Mass Function 


Ah 
(M ) 


Emerge 

(Gyr) 


J0755+4906 





06302 ± 


00213 


438 ±5 


-51 ±4 


5150.84949 ±0.00010 





550 


± 





025 


> 





81 


< 0.22 


J0818+3536 





18315 ±0 


02110 


170 ±5 


-201 ±4 


5151.90466 ±0.00030 





094 


± 





014 


> 





26 


< 8.89 


J0917+4638 1 





31642 ± 


00002 


145 ±2 


28 ±1 


3708.85755 ±0.00134 





102 


± 





003 


> 





27 


< 36.8 


J0923+3028 





04495 ± 


00049 


296 ±3 


2±2 


3818.67019 ±0.00002 





121 


± 





004 


> 





34 


< 0.13 


J1053+5200 2 





04256 ± 


00002 


264 ±2 


12 ±2 


3790.79731 ± 0.00004 





081 


± 





002 


> 





26 


< 0.16 


J1233+1602 





15090 ± 


00009 


336 ±4 


-35 ±3 


5268.84901 ±0.00014 





597 


± 





020 


> 





86 


< 2.14 


J1422+4352 





37930 ± 


01123 


176 ±6 


-195 ±8 


4596.92484 ±0.00068 





208 


± 





023 


> 





41 


< 42.9 


J1439+1002 





43741 ± 


00169 


174 ±2 


-55 ±1 


3882.68848 ±0.00079 





240 


± 





007 


> 





46 


< 54.7 


J1448+1342 3 








35 ±7 


-28 ±5 






















J1512+2615 





59999 ± 


02348 


115 ±4 


-41 ±4 


3879.48768 ±0.00119 





097 


± 





012 


> 





28 


< 171 


J1630+2712 





27646 ± 


00002 


218 ±5 


-118 ±3 


5009.72779 ±0.00016 





295 


± 





018 


> 





52 


< 15.6 


J2119-0018 





08677 ± 


00004 


383 ±4 


-28 ±4 


5008.88701 ± 0.00003 





501 


± 





016 


> 





75 


< 0.54 



References. — 1 iKilic et all (l2007bh : 2 iKilic et al.1 d2010bh 

Note. — 3 The measurements are formally consistent with no variation. 
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We present the best-fit orbital parameters in Table [2j Columns include orbital period 
(P), radial velocity semi-amplitude (K), systemic velocity (V^ yste mic), the time of spectro- 
scopic conjunction (the time when the object passes through km s _1 as it approaches 
the observer), mass function (see Eqn. 1 below), and minimum secondary mass (assuming 
% = 90°). The systemic velocities in Table |2] are not corrected for the WDs' gravitational 
redshifts, which should be subtracted from the observed velocities to find the true systemic 
velocities. This correction is approximately 3 km s _1 for our targets. 

We plot the best-fit orbits compared to the observed radial velocities in Figure [5J ex- 
cluding the two objects previously published. Follow-up observations were typically obtained 
over a pair of 3-day time baselines separated by one week, plus an original observation that 
extends the baseline up to 1488 days (not shown in the upper panel of Fig. E]). The ELM 
WD's short orbital periods combined with our long time baselines allow us to constrain the 
orbital periods accurately. 



4. DISCUSSION 
4.1. J0755+4906 

SDSS J075552. 40+490627. 9 exhibits the largest radial velocity variation in our sample 
of ELM WDs. The WD has a best-fit amplitude of 876 ± 10 km s , an orbital period of 
1.512 ± 0.051 hr, and a binary mass function of 



M 9 3 sin 3 i PK 3 
(MfTW = ^G = °- 550±0 -° 25Ms ' 



where % is the orbital inclination angle, M\ is the ELM WD mass, and M 2 is the com- 
panion mass. Although the inclination is unknown, we can use the distribution of possible 
inclinations to constrain the companion's mass. 

Given the observed orbital parameters, there is a 61% probability that J0755+4906's 
companion is a WD with <1.4 M and a 20% probability that the companion is a neutron 
star with 1.4-3.0 M . The remaining probability is for a companion mass >3 M . We 
estimate the most probable companion mass by assuming the mean inclination angle for a 
random stellar sample, i = 60°. For J0755+4906, the most likely companion is a 1.12 M Q 
object at an orbital separation of 0.7 R . Follow-up radio or X-ray observations are required 
to rule out a neutron star companion, but statistically the companion is most likely a massive 
WD. 



Short period binaries like J0755+4906 must eventually merge due to angular momentum 
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loss to gravitational wave radiation. The merger time is 



T 



(Mi + m 2 ; 

MiM 2 



1 1/3 



_p8/3 x 10 - 2 Gyr 



(2) 



where the masses are in M Q and the period P is in hours (ILandau fc Lifshitzl Il958h . For a 
companion mass of 1.1 M , J0755+4906 will merge in 170 Mvr. iKilic et all ( boiObl ) discuss 
the many possible stellar evolution paths for such a merging system. For J0755+4906, there 
is at least a 6% likelihood that the system contains a pair of WDs whose total mass exceeds 
the Chandrasekhar mass and thus will explode as a Type la supernova. On the other hand, 
the extreme mass ratio Mi/M 2 <0.21 means that mass transfer is likely to be stable and 
thus this binary will probably evolve into an AM CVn system. 



4.2. J0818+3536 

The best-fit orbital period for SDSS J081822.34+353618.9 is 4.396 ±0.51 hr, however the 
limited number of follow-up observations allow for aliases at 5.4 and 8.6 hr periods. Periods 
longer than 14 hr are ruled out. Thus J0818+3536 is a short period binary. Assuming the 
best-fit period, there is a 90% probability that the companion is a WD with <1.4 M and a 
5% probability that the companion is a neutron star with 1.4-3.0 M Q . For the most probable 
inclination angle i = 60°, the companion is a 0.33 M low mass WD at orbital separation of 



4.3. J0917+4638 



Previously published bv IKilic et all fl2007bf ). SDSS J091709.55+463821.8 has a binary 
orbital period of 7.5936 ± 0.0005 hr. Assuming an inclination angle of 60°, its companion is 
another low-mass WD with 0.35 M at an orbital separation of 1.6 R . 

Curiously, J0917+4638 and the two other ELM WDs with lower surface gravity in our 
sample (J1233+1602 and J2119— 0018) all show significant Ca II K absorption in their at- 
mospheres (Figure [2]). The extremely short timescale for gravitational settling of Ca in 
the WD photospheres means there must be an external source for the observed metals 
fcoester fc WilfeeiJliooa IKilic et al.ll2007af ). The ELM WDs are located far above the Galac- 
tic plane where accretion from the interstellar medium is unlikely. A plausible source of Ca 
is accretion from a circumbinary disk created during the mass loss phase of the giant. Since 
the ELM WDs all went through a common envelope phase with a companion, a left-over 
circumbinary disk is possible. 
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4.4. J0923+3028 



SDSS J092345.60+302805.0 is the brightest ELM WD in our sample, with g = 15.7 mag 
and a heliocentric distance estimate of 280 pc. J0923+3028 has a binary orbital period of 
1.079 ± 0.012 hr and a mass function of 0.121 ± 0.004 M . There is a 87% probability that 
the companion is a WD with <1.4 M Q and a 6% probability that the companion is a neutron 
star with 1.4-3.0 M Q . For the most probable inclination angle, i = 60°, the companion is 
a 0.44 M Q WD at orbital separation of 0.5 R^. Thus, J0923+ 3028 is the second shortest 
period double WD system after J1053+5200 (iKific et al.ll2010bl ). This system will merge in 
less than 130 Myr. 

J0923+3028 is the one ELM WD in our sample with a signifi cant proper motio n mea- 
surement, fi a cos5 = —4.2 ± 3.5 and fis = —25 ± 3.5 mas yr _1 (IMunn et al.l 120041 ) . The 
gravitational redshift of the WD is ~4 km s _1 , thus its true systemic veloci ty is — 2 ± 5 km 



The velocity components with respect to the local standard of rest (jSchonrich et al. 



20101) are U = 14 ± 6, V = -20 ± 8, W = -2 ± 6 km s _1 . J0923+3028 is clearly a disk star. 



4.5. J1053+5200 



Previously published bv iKilic et al.l fl2010bf ). SDSS J105353.89+520031.0 has a binary 
orbital period of 1.0224 ± 0.0005 hr. Assuming an inclination angle of 60°, the companion 
is another low mass WD with 0.33 M at an orbital separation of 0.4 R . This system will 
merge in less than 160 Myr. 



4.6. J1233+1602 

SDSS J123316. 20+160204. 6 is the lowest surface gravity WD known. Its surface grav- 
i ty, log q = 5.12 ± 0.07, i s comparable to a 2.5 M main sequence star with logg ~ 4.2 



( IGirardi et al.l I2002L 12004 ). However, J1233+1602's binary orbital parameters exclude the 



possibility of it being a main sequence star. 

J1233+1602 is a binary with a 3.6216 ± 0.0022 hr orbital period and a mass function 
of 0.597 ± 0.020 M . By assuming an edge-on orbit, % = 90°, we can place a lower limit 
on the companion mass. For a 2.5 M A star primary, the star must have at least a 2.44 
M companion at an orbital separation of 2.0 R . However, the radius of a 2.5 M main 
sequence star is 1.9 R . Thus the orbital separation of two A stars in this system is less 
than their summed radii. Even if the companion were a 2.44 M neutron star, the system 
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would still be in Roche-lobe overflow, for which we see no evidence. 

As a check, we combine the spectra near maximum blue-shift and red-shift into two 
composite spectra. If there is a contribution from the companion, it may be visible as 
second pair of lines or an asymmetry in the line profiles. We do not see any extra lines or 
asymmetries in the line profiles. We conclude J1233+1602 is a ~0.17 M© WD orbiting a 
much fainter companion. 

Given the observed mass function, there is a 58% probability that the companion is a 
WD with <1.4 M Q and a 21% probability that the companion is a neutron star with 1.4-3.0 
M Q . Assuming the most probable orbital inclination angle of 60°, J1233+1602's companion 
is a 1.2 M WD at an orbital separation of 1.3 R . Follow-up radio or X-ray observations 
are required to rule out a neutron star companion, but statistically the companion is most 
likely a massive WD. 

There is a 7% likelihood that this system contains a pair of WDs whose total mass 
exceeds the Chandrasekhar mass and thus will explode as a Type la supernova. The merger 
time for such a system is <2 Gyr. However, the extreme mass ratio Mi/M 2 < 0.20 suggests 
that mass transfer will be stable and this binary will likely evolve into an AM CVn system. 

4.7. J1422+4352 

SDSS J142200. 74+435253. 2 has the fewest follow-up observations yet its binary orbital 
period is reasonably well constrained. The best-fit period is 9. 10 ±0.27 hrs, with a significant 
alias at 15.3 hrs. Assuming the best-fit period, there is a 82% probability that the companion 
is a WD with <1.4 M and a 8% probability that the companion is a neutron star with 1.4- 
3.0 M . For the most probable inclination angle, i = 60°, the companion is a normal 0.55 
M WD at orbital separation of 2.0 R . 

4.8. J1439+1002 

SDSS J143948.40+100221.7 is a binary with a best-fit period of 10.498 ± 0.041 hrs. 
Despite having follow-up observations on 10 nights spaced over one year (see Fig. [5]) there 
remains a significant period alias at 18.576 ± 0.215 hr. Using the best-fit orbital solution, 
there is a 80% probability that the companion is a WD with <1.4 M and a 10% probability 
that the companion is a neutron star with 1.4-3.0 M . For the most probable inclination 
angle, i = 60°, the companion is a normal 0.62 M WD at orbital separation of 2.2 R . 
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4.9. J1448+1342 

SDSS J144801. 13+134232.8 is the one ELM WD for which we detect no significant 
velocity variation. The best-fit orbit (0.5 hr period) has a reduced y 2 identical to a constant- 
velocity fit, thus this system can have a velocity semi-amplitude no larger than K < 35 ± 11 
km s _1 . By comparison, the average semi-amplitude of the other eleven ELM WD binaries 
is 240 km s _1 . The required orbital inclination for J1448+1342 to have the average semi- 
amplitude of the other ELM WD binaries is i < 8°. 

J1448+1342 may be either a pole-on binary system or a single star. If the orbital 
inclinations of our non-kinematically selected ELM WDs are randomly distributed, we expect 
one of the twelve ELM WDs to have i < 23.5°; there is a 12% likelihood of finding one of 
the twelve systems with i < 8°. Thus it is possible that J1448+1342 is a pole-on binary 
system. The alternative is that J1448+1342 is a single star. It is intriguing that J1448+1342 
is the most massive object in our sample of 12 ELM WDs. Existing stellar evolution models, 
however, ca nnot explain 0.25 M WDs from single star evolution, even with extreme mass 



loss models (IKilic et al.ll2007d ). Additional observations are required to determine the nature 



of this ELM WD. 



4.10. J1512+2615 

SDSS J151225.70+261538.5 is a binary with a best-fit orbital period of 14.40 ± 0.56 
hr, and a significant period alias at 17.8 hrs. Assuming the best-fit period, there is an 
89% probability that the companion is a WD with <1.4 M and a 5% probability that 
the companion is a neutron star with 1.4-3.0 M . For the most probable inclination angle 
i = 60°, the companion is a 0.36 M low mass WD at orbital separation of 2.5 R . 



4.11. J1630+2712 

SDSS J163026.09+271226.5 has a best-fit orbital period of 6.6350 ± 0.0005 hr. The 
observations also allow for 5.16 and 9.17 hr periods. Assuming the best-fit period, there is 
a 77% probability that the companion is a WD with <1.4 M and a 11% probability that 
the companion is a neutron star with 1.4-3.0 M . For the most probable inclination angle 
i = 60°, the companion is a normal 0.70 M WD at an orbital separation of 1.7 R . 
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4.12. J2119-0018 

SDSS J211921. 96— 001825.8 is the second lowest surface gravity WD known, with log*? = 
5.36 ± 0.07. As with J1233+1602, J2119-0018's orbital parameters exclude the possibility 
of it being a main sequence A star. 

J2119— 0018 is a binary with a 2.0825 ± 0.0010 hr orbital period and a mass function 
of 0.501 ± 0.016 M . If the primary is 2.5 M main sequence star, its companion must be 
>2.2 M at an orbital separation of 1.4 R Q . This separation is smaller than the radius of 
main sequence stars. We do not see any evidence of a companion in our spectra. Thus 
J2119-0018 must be a -0.17 M Q WD. 

Given the observed orbital parameters, there is a 64% probability that the companion 
is a WD with <1.4 M Q and a 18% probability that the companion is a neutron star with 
1.4-3.0 M . Assuming the most probable orbital inclination angle of 60°, J2119— 0018's 
companion is a 1.04 M WD at an orbital separation of 0.9 R . Follow-up radio or X-ray 
observations are required to rule out a neutron star companion, but the companion is most 
likely a massive WD. 

This system will merge in less than 540 Myr. There is a 5% likelihood that the system 
contains a pair of WDs whose total mass exceeds the Chandrasekhar mass and thus will 
explode as a Type la supernova. However, the extreme mass ratio M\jMi < 0.23 suggests 
that mass transfer will be stable and this binary will likely evolve into an AM CVn system. 



5. CONCLUSION 



We present a complete sample of 12 ELM WDs with masses < 0.25 M in the HVS 
Survey. Eleven of the WDs are binary systems with orbital periods <14 hr; J1448+1342 
is the only WD for which we detect no velocity variation. Given the binary nature of the 
other WDs in our non-kinematically selected sample, J1448+1342 is quite possibly a pole-on 
system. Our observatio ns demonstrate th at the formation of <0.25 M ELM WDs requires 
close binary evolution (jMarsh et al.lll995l ). 



Two of our targets, J1233+1602 a nd J2119— 0018, h ave surface gravities lower than 
the previous record-holder J0917+4638 (IKilic et al.ll2007al ) and thus are the lowest surface 
gravity WDs currently known. Interestingly, all three of these low surface gravity ELM 
WDs exhibit strong Ca absorption in their spectra. The Ca is likely explained by on-going 
accretion from circumbinary disks. Other sources of accretion are much less likely, given the 
extremely short timescale for gravitational settling in the ELM WD photospheres and their 
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1-2 kpc distances above the disk plane. 



One possible source of circumbinary disk accretion is debris from former planetary 
systems (IDong et al. For the case of our ELM WD binaries, the systems must have had 

circumbinary planets that could survive the binary common envelope ev olution. Com mon 
envelope evolution may also lead to "second generation" planet formation (jPeretsll2010l ). the 
debris of which may be accreting onto the ELM WDs. In any case, the lowest surface gravity 
ELM WDs are interesting targets for future infrared observations. 

Six of our 12 ELM WDs (50%) will merge within a Hubble time (we exclude J1448+1342 
due to its unknown period); four of these mergers will happen in <600 Myr. In all cases the 
most likely binary companions are WDs, although neutron stars cannot be ruled out based 
on the available data. Two of the merger systems, J0818+3536 and J1053+5200, are likely 
to have unusually low-mass (<0.4 M p) WD com panions and thus are possible progenitors 
for single helium-enriched sdO stars (lHeberll2009l ). Four of the merger systems are likely to 
have C/O WD companions and thus will likely form extreme helium stars. However, three 
of the systems (J0755+4906, J1233+1602, and J2119— 0018) have such extreme mass ratios 
<0.23 that future mass transfer is likely to be stable and thus they will evolve into AM CVn 
systems. 

We expect that at least one of our ELM WDs has i > 85° and thus is an eclipsing 
system. ISteinfadt et al.l (120101 ) recently discovered the first eclipsing detached double WD; 
finding eclipsing WD binaries in our sample will provide important mass-radius constraints 
on He core ELM WD models. The most probable eclipsing candidates in our sample are 
the systems with the shortest merger times, namely J0923+3028, J1053+5200, J0755+4906, 
and J2119-0018. 

Our sample of 12 ELM WDs is unique in terms of its non-kine matic selection and 
its completeness. Pr evious discoveries o f ELM WDs in the SPSS (ILiebert et al 



Eisenstein et al.l 120061 ) and other surveys (IKawka et al.l 120061 ; iKawka fc Vennesl l2009f) have 



2004 



been u seful for studying the formation and future evolution of the individual systems (IKilic et al 
2010bl ). However, the spectroscopic selection biases present in the SDSS prohibit an accurate 
estimate of the space density of these objects. SDSS WDs were observed by different target- 
ing programs selected with different color, magnitude, and photometric qu ality cuts as well 



as di fferent sparse sampling rates that vary over the SDSS survey region (jEisenstein et al. 



20061 ). Our sample, on the other hand, is complete in magnitude, color, and spatial coverage 
and enables us to estimate the space density and merger rate of ELM WDs in the Milky 
Way. We explore possible links t o underluminous su pernovae, AM CVn stars, and other 
phenomena in a companion paper (IBrown et al.ll2010l ). 
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Fig. 5.— Data and best-fit orbits for the 10 new ELM WDs. Upper panels plot the 
heliocentric radial velocities vs. time. Lower panels plot the observations phased to the best- 
fit orbital solution (Table [2]), and also include the observation first used to identify the ELM 
WD candidate. The same vertical axis is used in all 10 plots. 
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Table 3. Radial Velocity Measurements 
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.2 


± 


20 


.7 


5277. 


.689765 


-287 


.2 


± 


31 


.7 


5277. 


,807903 


-155 


.8 


± 


24 


,1 


3818. 


,669260 


-22 


.3 


± 


13 


.4 


5204. 


.841296 


123 


.3 


± 


23 


.6 


5204. 


,849317 


364 


,1 


± 


29 


,1 


5204. 


.856366 


170 


.1 


± 


21 


.0 


5204. 


,864851 


-165 


.2 


± 


23 


.2 


5204. 


,871899 


-288 


.4 


± 


20 


.7 


5204. 


.880060 


-114 


.6 


± 


18 


.5 


5204. 


,887108 


101 


,8 


± 


24 


,9 


5204. 


,897757 


273 


.5 


± 


17 


.7 
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Table 3 — Continued 



Object 


HJD 






+2450000 


(km s _1 ) 



5204 


904794 


15.5 ± 18.6 


5204 


912596 


-198.1 ± 16.4 


5204 


927423 


-27.7 ± 16.8 


5204 


934460 


255.0 ± 11.6 


5204 


942365 


245.2 ±21.3 


5204 


949403 


49.5 ± 16.5 


5204 


957216 


-293.3 ± 21.2 


5204 


964265 


-248.4 ± 22.3 


5204 


972922 


2.3 ±20.3 


5204 


979960 


320.7 ± 15.0 


5204 


987969 


260.3 ±21.7 


5212 


845177 


171.1 ± 32.3 


5212 


858568 


183.6 ±45.1 


5221 


767105 


-248.2 ±28.5 


5221 


774790 


-307.1 ± 19.1 


5221 


790045 


153.6 ±29.9 


5221 


797082 


316.9 ±35.1 


5221 


804744 


85.6 ±26.1 


5221 


811793 


-170.3 ±21.6 


5221 


819339 


-318.7 ± 23.4 


5221 


826388 


-189.8 ± 35.6 


5221 


834351 


144.9 ± 17.4 


5221 


841400 


269.2 ±20.0 


5221 


848935 


133.3 ± 13.0 


5221 


855983 


-142.1 ± 15.5 


5221 


863657 


-282.4 ± 12.6 


5221 


870706 


-153.8 ± 18.7 


5221 


878275 


162.7 ±19.0 


5221 


885313 


300.7 ± 17.0 


5221 


892859 


129.4 ± 16.1 


5221 


899896 


-136.3 ± 17.6 


5221 


907431 


-266.4 ±20.3 


5221 


914468 


-177.7 ±29.9 


5221 


922038 


-15.5 ±22.4 


5221 


929086 


285.6 ± 18.0 


J1233+1602 5268 


822025 


-351.5 ± 10.6 
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Table 3 — Continued 



Object 


HJD 


Vhelio 




+2450000 


(km s -1 ) 



J1422+4352 



J1439+1002 



5268. 


.840648 


-157.5 ±13.1 


5268. 


.864607 


171.2 ±9.3 


5268. 


,887073 


301.6 ± 15.5 


5268 


.909573 


147.3 ±9.1 


5269. 


.793880 


290.8 ± 14.6 


5269 


,805316 


278.4 ± 17.4 


5269. 


.816670 


151.1 ± 15.8 


5270. 


.766544 


-344.9 ±14.1 


5270. 


,775294 


-347.4 ± 20.0 


5276. 


,854668 


97.5 ± 16.4 


5276. 


.866474 


191.1 ± 12.1 


5276. 


,875895 


298.3 ± 12.6 


4596. 


,887690 


-290.1 ± 10.9 


5268 


,969670 


-364.8 ± 12.5 


5268. 


.978929 


-331.1 ± 9.8 


5268 


.988397 


-352.3 ± 13.6 


5270. 


.854102 


-373.6 ± 14.8 


5384. 


,771754 


-93.0 ± 15.4 


5384 


,793095 


-36.8 ± 15.6 


5385. 


,726490 


-314.3 ±21.4 


3882 


885571 


—8.6 ± 7.9 


4922. 


,892856 


-21.0 ±6.9 


4922. 


.892856 


82.0 ±4.5 


4922. 


,977535 


116.4 ±5.6 


4923. 


.020372 


137.9 ± 18.8 


4923. 


,973820 


-126.2 ±7.4 


4924. 


,020430 


-191. 7± 18.9 


5267. 


.943738 


-56.8 ±14.4 


5267. 


.969654 


-50.0 ± 11.2 


5268. 


.036672 


-15.3 ±25.6 


5268. 


,957759 


111.4 ±7.0 


5268. 


,962158 


110.6 ±5.1 


5269. 


,971779 


-147.0 ± 14.4 


5269. 


,977266 


-165.6 ±11.1 


5270. 


.959442 


-189.1 ±5.7 


5270 


,979860 


-178.5 ±9.1 
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Table 3 — Continued 



Object 


HJD 


Vhelio 




+2450000 


(km s -1 ) 



J1448+1342 



J1512+2615 



5271. 


.011239 


— 119.. 


3 ±9.1 


5275. 


.934512 


120, 


.3 


± 


10.2 


5277. 


,886834 


-216, 


.8 


± 


10.5 


5384. 


.660342 


-246 


.8 


± 


10.3 


5385. 


.755771 


202 


.7 


± 


47.3 


4237. 


,919202 


-21 


.6 


± 


15.6 


5276. 


.885957 


-38 


.5 


± 


19.0 


5276. 


.892265 


-0, 


.1 


± 


25.0 


5276. 


.900298 


-22 


.9 


± 


17.6 


5276. 


,906641 


-39 


,8 


± 


23.6 


5277. 


.023591 


-182 


.0 


± 


69.6 


5277. 


,896227 


-3 


.2 


± 


25.8 


5277. 


,902686 


-14 


.6 


± 


27.1 


5277. 


,911552 


-49 


.7 


± 


32.8 


5277. 


,917988 


4, 


.5 


± 


20.8 


5277. 


.990376 


-53 


,4 


± 


22.3 


5277. 


,996811 


-82 


.0 


± 


27.8 


5278. 


.004983 


-35 


.8 


± 


20.4 


5278. 


.011256 


-23 


.8 


± 


16.6 


5384. 


.666893 


-0, 


.3 


± 


16.7 


5384. 


752430 


-56 


,4 


± 


20.0 


5385. 


,762534 


-12 


.3 


± 


32.9 


3879. 


.873052 


-127 


.7 


± 


11.2 


5267. 


,949987 


56, 


.7 


± 


17.4 


5267. 


.974827 


87. 


.6 


± 


19.0 


5267. 


,997247 


49 


,5 


± 


18.6 


5269. 


,017902 


-96 


.7 


± 


14.4 


5269. 


.023840 


-101 


,2 


± 


17.1 


5270. 


.006995 


-74 


.1 


± 


19.5 


5271. 


.002158 


66. 


.2 


± 


16.6 


5271. 


.017749 


1 


.0 


± 


28.4 


5271. 


.024775 


53 


.1 


± 


26.8 


5276. 


.915247 


42. 


.9 


± 


17.9 


5277. 


.013574 


129 


.0 


± 


31.5 


5277. 


.943822 


-154 


.4 


± 


25.1 


5384. 


,675640 


-189 


.8 


± 


22.7 



Table 3 — Continued 



Object 


HJD 


Vhelio 




+2450000 


(km s -1 ) 



J1630+2712 



J2119-0018 



5384.763249 
5385.747861 
5009.892692 
5267.958434 
5267.983829 
5268.011944 
5269.005570 
5269.985526 
5270.988168 
5276.925057 
5276.999529 
5384.809434 
5384.868187 
5385.781067 
5385.842754 
5008.884992 
5152.578978 
5152.591870 
5153.581786 
5153.590153 
5153.597282 
5153.605372 
5153.612500 
5153.628726 
5153.636966 
5153.645148 
5153.654661 
5153.663884 
5157.626210 
5157.635329 
5157.644287 
5157.653291 
5157.662260 
5157.672131 
5158.586060 
5158.594878 



-127. 


.7 


± 


15.2 


22. 


.7 


± 


28.1 


-239 


.7 


± 


18.9 


-47. 


.1 


± 


18.0 


24. 


.0 


± 


14.3 


103 


.0 


± 


14.5 


-280 


.4 


± 


13.7 


17. 


.1 


± 


15.1 


-80 


.9 


± 


14.5 


-133 


.0 


± 


19.3 


-332 


.8 


± 


29.4 


-345 


.1 


± 


15.2 


-197 


.8 


± 


31.1 


126 


.8 


± 


13.6 


-58 


.9 


± 


34.1 


-15 


.2 


± 


18.7 


7. 


.8 


± 


10.5 


266 


.2 


± 


10.8 


-203 


.5 


± 


16.9 


-319 


.1 


± 


21.7 


—396 




± 


15.3 


-329 


.3 


± 


18.1 


-159 


.7 


± 


15.1 


306 


.2 


± 


23.9 


357 


.4 


± 


17.9 


312 


.6 


± 


18.3 


154 


.2 


± 


28.5 


-115 


.9 


± 


29.2 


383 


.6 


± 


31.6 


343 


.2 


± 


16.8 


265 


.7 


± 


17.8 


26. 


.0 


± 


28.6 


-139 


.8 


± 


25.8 


-427 


.6 


± 


17.2 


337 


.8 


± 


21.8 


312 


.5 


± 


16.8 
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Table 3 — Continued 



Object HJD 


^helio 


+2450000 


(km s- 1 ) 


5158.603454 


150.1 ±30.5 


5158.613036 


-173.7 ±25.2 


5158.625477 


-383.8 ± 17.5 


5158.636229 


-443.9 ± 18.3 



